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£3 Abstract 

In 1995 Magnus [15] posed a conjecture about the asymptotics of the recurrence coefficients 
of orthogonal polynomials with respect to the weights on [—1, 1] of the form 

(l-zni + ^lxo-zPx/ 5 ' fo^t" 1 ^)' 

[^4, for x G [xo, 1] , 

with A, B > 0, a, (3, 7 > —1, and xq G (—1,1). We show rigorously that Magnus' conjecture 
ON is correct even in a more general situation, when the weight above has an extra factor, which 

is analytic in a neighborhood of [—1, 1] and positive on the interval. The proof is based on 
the steepest descendent method of Deift and Zhou applied to the non-commutative Riemann- 
Hilbert problem characterizing the orthogonal polynomials. A feature of this situation is that 
the local analysis at xq has to be carried out in terms of confluent hypergeometric functions. 
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1. Introduction and statement of results 



1.1. Introduction 

A. Magnus considered in [15! a weight function which is smooth and positive on the whole 
interval of orthogonality up to a finite number of points where algebraic singularities occur. 
His primary goal was to investigate the influence of these singular points on the asymptotic 
behavior of the recurrence coefficients of the corresponding orthogonal polynomials (general- 
ized Jacobi polynomials). Based on numerical evidence, he conjectured explicit formulas for 
the asymptotics of these coefficients (as the degree of the polynomial grows) for the weights 
of the form 

[A, for x e [x , 1J , 

with A, B > and 7 > —1, and xq € (—1, 1). This weight combines at a single point 
both an algebraic singularity and a jump. 

So far, Magnus' conjecture has been confirmed rigorously in some special cases (see below); 
our main goal is to establish it in its full generality, and even to extend it further. Namely, 
we consider polynomials that are orthogonal on a finite interval [—1,1] with respect to a 
modified Jacobi weight of the form 



c,7 



(x) = (l-x) a (l + x)P\xo-xf<h(x)Z c (x), a; €[-1,1], (2) 



where xo € (—1, 1), a,/?, 7 > —1, h is real analytic and strictly positive on [—1, 1], and E c is 
a step-like function, equal to 1 on [—1, xo) and c 2 > on [xo, 1]. 

The proof is based on the nonlinear steepest descent analysis of Deift and Zhou, introduced 
in [7] and further developed in [SJ [61 [9] , which is based on the Riemann-Hilbert characteri- 
zation of orthogonal polynomials due to Fokas, Its, and Kitaev |10| . A crucial contribution 
to this approach is [13], where the complete asymptotic expansion for the orthogonal poly- 
nomials with respect to the Jacobi weight modified by a real analytic and strictly positive 
function was obtained (in notation (JlJ), A = B and 7 = 0). The first application of this 
technique to weights with a jump discontinuity is due to [13], where the authors considered 
an exponential weight on M with a jump at the origin. 

Let P n (x) = P n (x; ui Cj7 ) be the monic polynomial of degree n orthogonal with respect to 
the weight w cn on [—1,1], 

1 

P n (x; w cn )x k w cn (x) dx = 0, for k = 0, 1, . . . , n — 1. 

1 

It is well known (see e.g. [16 ) that {P n } satisfy the three-term recurrence relation 

P n+ l(z) = {z- b n )P n (z) - a 2 n P n ^(z). (3) 
The central result of this paper is: 
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Theorem 1. The recurrence coefficients a n and b n of orthogonal polynomials corresponding 
to the generalized Jacobi weight Q have a complete asymptotic expansion of the form 



1 _ A k {n) 

2 ^ n k 

k=l 



£ 

k=l 



B k {n) 



ir 



as n —* oo, where for every k € N the coefficients A k {ri) and Bj.(n) are bounded in n. In 
particular, 



Ai(n) 
Bi(n) 

where 



' 7 2 log 2 c 
— H ^— cos 

4 7T 2 



2n arccos xo + 2 



log c 



7T 



log ( 4n\/l - x 2 



e 



1 — X 



2. IT , lo g c 



+ 



COS 



7T- 



log c . 



(2n + 1) arccos xq + 2 log ( 4ny 1 — x 



e 



(4) 



(5) 



e 



a + — ) 7T— (a + /? + 7) arccos xq — 2 arg T 



7 . log c 



7T 



7 .logc\ y 7 ! ~ x o j} ^°gh(t) dt 



-arg[ 

Z 7T 



7T 



(6) 



-1 VT^t-xo , 

and j denotes the integral understood in terms of its principal value. 
Remark 1 . We can rewrite the result of this theorem as 



1 M 

2 ~ n 

2M 



cos 



2n arccos xq — 2fi log I AnJ 1 — Xq ) — 



+ 



1 



cos 



n 



as n — > 00, where 



(2ra + 1) arccos xo — 2/i log ( 4n^/ 1 — Xq 
log c 



M 



e 



+ 



n- 



+ IJ- 2 , 



vr 2 V 4 

and defined by (j6J). 

A comparison of these formulas with those in |15j (setting h(x) = .B and c 2 = A/B) shows 
that Magnus' conjecture on the asymptotic behavior of the recurrence coefficients holds true 
for weights of the form ([T]). Observe that this is a slight extension of the original statement 
of Magnus: (i) we allow for an extra real analytic and strictly positive factor h in the weight, 
and (ii) we can replace the error term o(l/n) in [15] by a more precise C(l/n 2 ). 
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Theorem [T] generalizes some previous known results about the asymptotics of the recur- 
rence coefficients. To mention a few, weight wi t o was considered in [T5], w\ n is a particular 
case of the weight studied in |17] , and w Ct o was matter of attention in [11] . 

The proof is based on the steepest descent analysis of the non-commutative Riemann- 
Hilbert problem characterizing the orthogonal polynomials P n . In theory, this approach 
allows to compute all coefficients and in Q-Q. However, the computations are 
cumbersome and their complexity increases with k. 

Most of the steps of the steepest descent analysis below are standard and well described in 
the literature, see e.g. (5j El H3] • The main feature of the situation treated here in comparison 
with the classical Jacobi weight is the singularity of the weight of orthogonality at Xq. The 
local behavior of P n 's at xq is described in terms of the confluent hypergeometric functions. 
Such functions appeared already in the Riemann-Hilbert analysis for the weight w c q in |11] 
and [13], and will work also in the general situation studied here. However, the parameter 
describing the family of these functions is complex; its real part depends on the degree of 
the algebraic singularity 7 and its imaginary part is a function of the size of the jump c 2 in 
the weight w cn . Also, for c = 1 these confluent hypergeometric functions degenerate into the 
Bessel functions, in correspondence with [14 . 

Interestingly enough, the confluent hypergeometric functions appear in the scaling limit 
(as the number of particles goes to infinity) of the correlation functions of the pseudo-Jacobi 
ensemble, see [3J. This ensemble corresponds to a sequence of weights of the form 

(1 + x 2)-n-Rc(s) e 2Im( S )arg(l+ J z) ) % £ R? ^ 

where n is the degree of the polynomial and s a complex parameter. The connection between 
both problems becomes apparent if we perform the inversion x 1— ► 1/x in dTI); this creates at 
the origin an algebraic singularity with the exponent Re(s) and a jump depending on Im(s). 

In the next Section we state the Riemann-Hilbert problem for the orthogonal polynomials 
and perform the steepest descent analysis; as a result, Theorem [T] is proved in Section 3. For 
the sake of brevity, the description of the standard steps is rather sketchy; an interested reader 
may find the missing details in the literature cited above, and especially in [TT]. However, 
the local parametrix (the Riemann-Hilbert problem in a neighborhood of the singularity) at 
xq is analyzed in full detail. The same problem has appeared very recently in an independent 
work of Deift, Its and Krasovsky [4] on the asymptotics of Toeplitz and Hankel determinants. 

2. The steepest descent analysis 

2.1. The Riemann-Hilbert problem and first transformations 

Following Fokas, Its and Kitaev [10 , we characterize the orthogonal polynomials in terms 
of the unique solution Y of the following 2x2 matrix valued Riemann-Hilbert (RH) problem: 
for n £ N, 
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(Yl) Y is analytic in C \ [-1, 1]. 

(Y2) On (— 1, xq) U (xo, 1), Y possesses continuous boundary values Y + (from the upper half 
plane) and Y_ (from the lower half plane) , and 



Y + (x) = Y_(x) 



1 W Cl7 (x) 

1 



(Y3) Y(z) = (I + C(l/z)) z nas , as z — » oo, where all terms are 2x2 matrices, I is the 



identity and 03 



1 



is the third Pauli matrix 



(Y4) if (C, s) e {(-1,(3), (x ,7), (1, «)} then for z -> C, 2 6 C\ [-1, 1], 



Y(z) = 



1 I* -en 

/ I- -CIV' 

1 io g |z-ch 

1 log^-ClJ 

1 ') 

1 1 ' 



if s < 0; 
, if s = 0; 
if 5 > 0. 



Standard arguments (see e.g. p3|) show that this RH problem has a unique solution given by 



Y(z,n) 



Pn (z) C (P n W c>7 ) (z) 

-2mk 2 n _ x P n _ x (z) -2mkl_ x C (P n _iio C)7 ) (z) 



(8) 



where P n is the monic orthogonal polynomial of degree n with respect to w crn p n {x) = 
Pn(x;w c ,-y) is the corresponding orthonormal polynomial, 

Pn(x) = k n P n (x), 

where k n > is the leading coefficient of p n , and C (•) is the Cauchy transform on [—1, 1] 
defined by 

1 f l fix) 
2ni J_i x — z 

Note that Y and other matrices introduced hereafter depend on n; however, to simplify 
notation we omit the explicit reference to n. 



C(f)(z) 



dx . 



1 In what follows, for a G C \ {0} and b € C, a btT3 we use the notation 

a" 
.0 l/a b 



dcf ( (X 

a 1 = 
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The first transformations of the Deift-Zhou steepest descendent method are standard, 
and up to slight variations match those described in subsection 2.2]. Hence, we will omit 
the details, highlighting basically the differences with the cited reference. The reader should 
keep in mind also that in [TT] the analysis is made for a singularity fixed at xo = 0; however, 
extending it to the general case of xo £ (—1,1) is straightforward. 

We start by defining 



T(z) 



Y(z)<p(z) 



0) 



where 



is the conformal map from C \ 

V. 



lf (z) = z+v / z 2 -l (10) 
-1,1] onto the exterior of the unit circle, with the branch of 



z* — 1 analytic in C \ [—1,1] and that behaves like z as z — > oo. 
In order to perform the second transformation we need to extend the definition of the 
weight of orthogonality to a neighborhood of the interval [—1,1]. By assumptions, h is a 
holomorphic function in a neighborhood U of [—1,1], and positive on this interval. For any 
Jordan arc T, intersecting [—1,1] transversally at xq and dividing U into two connected 
components, we denote by S5 its intersection with the upper half plane, and by T,q its 
intersection with the lower half plane, oriented as shown in Figure [l^ Contours S5 U T,q U R 
divide U into four open domains ("quadrants"), that we denote by Q± R as depicted. Finally, 
let Q L (resp., Q R ) be the connected component of U \ T containing —1 (resp., +1). 



Q 



R 



Q 1 - 



Q 1 - 



Figure 1: Division of the neighborhood of [—1, 1] in four regions. 



Now we set 



w« s » w (i - (i + ,y> x h - % • * e QL R \\;™' - 1 ! • 

\(z-xo) 7 , z € Q K \ [1, +00) , 



111) 



where the principal branches of the power functions are taken. In this way, w is defined and 

"dcf 



holomorphic in U = U\ ((—00, —1] U [1, +00) U V), and w(x) > for x E ( — 1, 1) \xq. Setting 

r l zeQ L 
c 2 z G Q R , 
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we extend also 

w c ,-y(z) = w{z) E c (z) , (12) 

to a holomorphic function in U. 

We describe now the next transformation consisting in opening of lenses or contour de- 
formation. 




Figure 2: First lens opening. 



We build the four contours Sj lying in U (except for their end points) such that Si and 
£3 are in the upper half plane, and Si and S2 are in the left half plane, and oriented "from 
— 1 to 1" but now through xo (see Fig. [2]). This construction defines three domains: the 
inner upper domain, bounded by [—1, 1] and the curves Si and S3; the inner lower domain, 
bounded by [—1,1] and the curves S2 and S4, and finally the outer domain, bounded by 



curves Sj and containing the infinity. Denote S 
Using the matrix T from d9j) we define 



[-MluUjUiE 



'Jfe- 



S(z) 



T(z) 
T(z) 

T(z) 



for z in the outer domain, 

1 (A 

1 -2 / \ 1 I ' f° r z m ^ ne mner upper domain, 



for z in the inner lower domain. 




(13) 



Then S is a solution of a new RH problem, now with jumps on S, that are easy to compute 
explicitly. The uniqueness is guaranteed if we impose the additional local requirement: as 

z — ► xq, z £ C\S, 



for -1< 7 < 0, S (z) 



O 



1 \z — Xq\ 
1 \z — Xq\ 



7 \ ' 



as z — > xq, 



7 



for 7 = 0: 



S(z) = { 



for 7 > 0: 



O 
O 



1 log|z-xo|\ 

l iog|z-xo|y ' 

log I Z — Xq I log I z — Xq I 
log|z — xq| loglz — Xq\ 



as z — » xq from the outer domain, 
as z — » xq from the inner domains, 



S(z) = 








1 1 
1 1 



-7 



|Z — Xq\ ' 1 
Iz-XoT 7 1, 



as z — > xq, from the outer domain, 
as z — » xq, from the inner domain. 



Outer parametrix 

In the next step, which is also standard, we build the so-called outer parametrix for the 
RH problem for S in terms of the Szego function D(-, u> Cj7 ) corresponding to the weight u> Cj7 . 
Namely, we construct the 2x2 matrix N that satisfies 



(Nl) N is analytic in C\ [-1, 1]; 



-w cl (x)" 



(N2) N+(x) = N_(x) 
(N3) N(z) = I + O (1/z), as z -> oo. 
The solution is given by 







x G (-l,x ) U (x , 1); 



N(z)^D%A(z)D(z,w cn )-° 3 , 
and we will describe the three factor appearing in the r.h.s. of ( |14| ). Matrix A is 

<p(z) 1 



(14) 



A(z) 



A n A 12 
-A 12 An 

A 12 (z) 



A u (z) 



V2(z 2 -1 



1/4 ' 



-1/2 



(15) 



V2(*2_i)V4 ^(z) 



with the main branches of the roots, in such a way that An is analytic in C\ [—1, 1] with 
An (z) — > 1, and A12 (z) — ► 0, as 2; — > 00. The Szego function D(-,w ca ) for uj Cj7 is given by 



D(z, w Cjl ) = D(z, h)D(z, wi n )D(z, S c ) , 



(16) 
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where 
D(z, h) = exp 

and 
with 



(z-l) Q/2 (* + !)"/*(* -so) 



/3/2 



(/9 (a+/3+ 7 )/2( z ) 



-D(z,H c ) = c exp —A log 



1 — zxq — iy z 2 — 1 



z - x 



def .logc 
A = z . 

7T 



^7/2 
(17) 

(18) 
(19) 



where we take the main branches of (z — l) a ^ 2 , (z + l)^ 2 , (z — xo) 7 ^ 2 and \/ z 2 — 1 that are 
positive for z > 1, as well as the main branch of the logarithm (see |11| section 2.3] for a 
detailed computation). Finally, 



Doo = D(oo, w c 



r c D (OO, h) 2~( Q +' 3 +T)/ 2 g^arcsinxo > n 



i ~ / \ • I \ ■ . M / _ ( v. (^^) 

Some of the properties of this function are summarized in the following lemma: 



Lemma 2. The Szego function D(-, w) for the weight w defined in (11) exhibits the following 
boundary behavior: 

lim D(z,w) = \/w(x) exp ( ±i$(x) ) , 



(21) 



with 



$ (x) 
$ (x) 



z-»a:e(-l,l), 

±lmz>0 



vra a + /? + 7 Vl - x 2 f 1 log h(t) dt 

arccos x -r , , 

2 2 2ir J_ iy /i—pt-x' 



$(x) + -1 < X < x 
$ (x) , x < x < 1 



(22) 
(23) 



Furthermore, for the step function H c , 



lim 

±lmz>0 



D(z, S c ) = v^M exp p. 



log c 



7T 



log 



1 - X X + \J{1- X 2 ) (l - Xq) 



X — Xq 



and 



D(z,E c ) = c ±n l0g ^~^J [1 + o(l)] , dsz^xo, ±lmz>0. 



(24) 
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The proof of this lemma is similar to Lemma 7] , up to the difference that our jump 



here takes place at a generic point xq, and that w{z) (see (11)) has an extra factor which 
makes w{z) non-analytic across S5 U T,q. 

The main purpose for constructing N is that it solves the "stripped" RH problem, ob- 
tained from the RH problem for S by ignoring all jumps asymptotically close to identity. 
Unfortunately, this property of N is not uniform on the whole plane: the jumps of N and S 
are no longer close in the neighborhoods of ±1 and xq. The analysis at these points requires a 
separate treatment, called local analysis, that we perform next. The outline of this analysis at 
each point is the following: take a small disc centered at the point and build a matrix-valued 
function (local parametrix) that: 

(i) matches exactly the jumps of S within the disc, and 

(ii) coincides with N on the boundary of the disc, at least to an order o(l), n — > 00. 
2.3. Local parametrix 

We fix a 5 > small enough such that discs = {z G C : \x — Q\ < 5} , ( G {— 1, xo, 1} 
are mutually disjoint and lie in the domain of analyticity of the function h. We skip the 
details of construction of the local parametrices P-ti at z = ±1 and refer the reader to |14] . 

For the local parametrix at the jump we need to build a 2 x 2 matrix-valued function 
P xo = Po in U X() \ T, that satisfies the following conditions: 

(Pol) Po is holomorphic in U XQ \E and continuous up to the boundary. 
(Pq2) Pq satisfies the following jump relations: 



1 

w cn (z 



for z G U X0 n ( |J Si I \ {x }; 



P +(x) = P -(x) i_ 'J I' for Z£ U xo n((-l,x )U(x Q ,l)) 

(P 3) P (z)N- 1 (z) =1 + (1/n), as n -» 00, uniformly for z G dU XQ \H. 
(Po4) Po has the same behavior than S as z — » xq, z G U Xo \H. 

Following a standard procedure, we obtain the solution of this problem in two steps, 
getting first a matrix P^ that satisfies conditions (Pol, Po2, Po4). After that, using an 
additional freedom in the construction, we take care of the matching condition (Po3). 

We define an auxiliary function W, holomorphic in U \ R. In the next formula we under- 

stand by I h (z) (1 — z) a (1 + z) (z — xo) 7 c I the holomorphic branch of this function in U\ 
((—00, xq] U [1, +00)), positive on (xq, 1). Analogously, ( h (z) (1 — z) a (1 + zy (xq — z) 1 c ) 
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stands for the holomorphic branch in U \ ((— oo, —1] U [xq, +oo)), positive on (—1, xq). With 
this convention we set 



W{z) 



(h{z){l-z) a {l + zf{z-x Q y C ) 

(h{ z ){i- z) a {i + zf ( XQ - z y ^ 



1/2 



1/2 



z G Q+UQ L , 
z G U Q?. 



It is easy to see from (11 >— ( 12 ) that 



W 2 (z) 



w cn {z)e~^ i c~\ 
w cn (z) e 7? "c, 
w cn (z) e~ 77 "c, 
w cn (z) e 77 ™c _1 , 



z G Q\, 
z G Q R . 



This shows that W satisfies the following jump relations: 

W+{x) 



W- (x) e^ n , 



— 1 < X < Xo, 
Xq < X < 1, 



and 

Moreover, 



W+ (z) = e^ n/2 (z), z G S 5 U S e 



W + (x) 



y/w cn (x) ce* 7 2 = ^/lui^, (x) 



- ,- 7T 

ce 2 . 



and 



(x) W- (x) 



ce 



— 1 < X < Xq, 
Xq < X < 1, 



J u; Cj7 (x) c, — 1 < x < xo, 
1 w Cj7 (x) c _1 , xo < x < 1. 
We construct the matrix function Pq in the following form: 



'0 U 



E n (z)pW (z)W(z)- <73 9 9(^ 



(25) 



(26) 



(27) 



(28) 



(29) 



where E n is an analytic matrix-valued function in U xo (to be determined). Matrix pW is 
analytic in U XQ \ E. Denote by 



Ji 







Ja 



l/c OJ' 

1 

e 7 ™c 1 



1 

e -7Tt c -l 1 

l/c 
-c 



J 5 



J 3 = «/ 7 



-,1771-/2 







e"^/ 2 

1 

e^c 1 



(30) 
(31) 
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Using the properties of W and (p it is easy to show that 



PL 1} (x) 



and 



Pi 1} M 



(z) <( 



J 5 , x G (a?o — #,x ), 
yi, x G {x ,x + 5), 

J 4 , z g Si n C4 \ {x }, 

J 6 , z G S2 n u xo \ {x }, 

J 2 , Z G S3 D f7a; \ {x }, 

J 8 , Z G S 4 Pi f7a; \ {x }, 



(32) 



(33) 



and, as W has a jump on S5 U Sg, by (26), we have two additional jumps on S5 U Se 

P? (*) 



P« (z) 



j 3 , z g s 5 n u XQ \ {x }, 

J 7 , z G s 6 n U X0 \ {x }. 



(34) 



Taking into account that W (z) = O (\z — xo| 7//2 ) and <p(z) = O (1) as z 
conclude also from (Po4) that P^ 1 ) has the following behavior at xq\ as z - 
C\(SUS 5 US 6 ), 



for 7 < 0: 

for 7 = 

P«(z) 

for 7 > : 
P«(z) = 



PW(; 



O 



\z - x | 7/2 

\z - Xq\ i/2 



O 
O 



o 
o 



log|z-x | log|z-x | 
log|z-x | log|z-x | 
1 log|z-x | 
1 log|z-x | 



\z — Xq[ 1 ^ 2 \z — Xq\ 7//2 ' 

\Z - X \ l/2 \Z - X \~ l/2 j 

\z — xo|~ 7 ^ 2 \z — xop 7 ^ 2 

I z — Xq I ~ I z — Xq I ~ 7 ^ 2 



\z - x | 7/2 
\z - x | 7/2 



from inside the lens, 
from outside the lens, 



-> xo, we 

X , Z G 



(35) 



(36) 



from outside the lens, 
from inside the lens. 



(37) 
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Figure 3: Auxiliary contours T. 

In order to construct P^ 1 ^ we solve first an auxiliary RH problem on a set T = Uj=i Tj 
of unbounded oriented straight lines converging at the origin, like in Fig. [3] More precisely, 

T 1 = \te™l 2 : t > 0} , T 2 = |te 3i7r / 4 : t > o} , T 3 = {-t : t > 0} , T 4 = {te 5i ^ 4 : t > o} , 
T 5 = {te 3in/2 : t > 0} , T 6 = |te-^* /4 : t > o} , T 7 = {t : t > 0} , T 8 = {e int/4 : t > o} . 

These lines split the plane into 8 sectors, enumerated anti-clockwise from ® to ® as in Fig. [3] 
We look for a 2 x 2 matrix valued function \I/ (z), satisfying the following conditions: 

(<3>1) * is analytic in C\I\ 

(^2) for k = 1, . . . , 8, \& satisfies the jump relation = on T^, with given 



by (30) and (31). 



(vl/3) the behavior of \& as £ — > is obtained from that of P*- 1 -* at x$ by replacing (z — xq) 
with £. Now the region "inside lens" correspond to ®U®U©U® and the region 
"outside lens" corresponds to ® U ® U © U ®. 

We construct \l/ explicitly using the confluent hypergeometric functions 

4> (a, 7 + 1; C) = 1F1 (a; 7 + 1; C) and V (a, 7 + 1; = z~\F Q (a, a - 7; -; -I/O , 

that are solutions of the confluent hypergeometric equation Qw" + (7 + 1 — C) w' — aw = 0, 
see [U formula (13.1.1)]. Namely, let 

G (a, 7; C) = C /2 <P (a, 7 + 1; H (a, 7; C) = C 7/ V (a, 7 + 1; e^ /2 , (38) 
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they form a basis of solutions of the confluent equation (see e.g. [TJ formula (13.1.35)]) 



4C V + 4(w' + [- 7 2 + 2C (7 + 1 - 2a) - C 2 



w 



0. 



(39) 



We can relate G and H with the Whittaker functions: G(a, 7; z) = M Kfl (z) / \fz and H (a, 7; 2) = 
W K|M (^)/v^ with [i = 7/2 and « = 1/2 + \i - a (see formula (13.1.32)]). 

In general, G (a, 7; Q and -ff (a, 7; () from (39 ) are multi- valued, and we take its principal 
branch in — | < arg (£) < For these values of ( we define 




-H(\ + lr,() 
%^G(l + A + i,7;C) y^H(l + \+lr,C) 



e 4 



<J3 



*(C) 



*(C)«, 


for £ G ® 


* (C) -Wl-^, 


for C G © 


*(C)^JlJ 2 J3, 


for c e ® 


* (C) ^8 Jl Jl h J±\ 


for c e © 


* (C) Jf 1 ^, 


for C G © 


*(C)^\ 


for C G © 


*(0, 


for C G © 


* (0 J 8 , 


for C G © 



then \l/ has the jumps across T specified in (^2)- Explicitly, \I/ (£) 



_r(i+A+i 



. £(iZ^±l) rr /j _ A , 7 />-^ 



-^(1 + A + 2, 7 ;C) 



£ feSf ) G(A + i,7;C)e- 



7 7T i 

~2~ 



r(i-A+; 



£te+il G (1 + a + 2 



r(7+i) 
r(i-A+ 



2 ' 



T' 
e 2 



T(Fx p(l-A + i,7;Ce-) 
F(2-A, 7 ;Ce- m ) 



T^fG(A + i,7;C) 



%gflG(l + A + i,7;C) 



r(i+A) " ^ " 1 2' 
F(2-A, 7 ;Ce- m )e 



7; Ce J e 2 



V 



r(i+A) ^ V 1 A + 2 '7^e 
^(2-A, 7 ;Ce- m ) 



(40) 



(41) 
(42) 
(43) 



e 4 



<T3 



(44) 
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. E^±|i H (i-x + ir, C e-) -H(x + l r , c) 



H{1 



A, 7; Ce 7 ™) e" A - r r(i- + A) ) g (l + A + 2 , 7; () 



77ri 



(45) 

«^(l + A +il7 ;0 ^ H{l + X + lrX) )* (46) 

^(i + a + ^;C) «, (1 + A + i,7;C)j e43 ' CG@ ' (47) 



-c" ir |#ff (1 - A + lr, Ce-«) (A + |,7; C) 



1 H (1 - A, 7; Ce-") £ ^f F (1 + A + i, 7; C) 



e— CT3 , C G ®- 



(48) 



Proposition 3. T/ie solution of the RH problem (^l), ( l I / 2), (^S) is given by (40) and 
det*(z) = 1, for z eC\T. 

Proof. If we take the branch cut across arg£ = —tt/2 oriented towards the origin (we 
consider —tt/2 < arg£ < 3ir/4), we have that the matrix \l/ has on this cut the following 
jump (using (19)): 

*+(C) = *-(C)J 5 , Cer 5 , (49) 
*+(0 = *-(C)( e -*^ J' Cerv (50) 

Using the following relations (see [E2 appendix: formulas (7.18), (7.30), (7.27)]), 

4> (a, b; e ±27Ti z) = <f> (a, b; z) , (51) 

V> (a, 6; e 2 "z) = e" 2 -> («, &; *) + e """ Y (a) T (1 + a - ft/ (& " °' 1; ^ ^ 
^ (6 - a, 6; e™z) e z = e^V (a, b; z) + ^je~ i7r ^<f> (a, 6; Z ) , 

r( s )r(i- s ) 



e vna _ g- 

and, combining the last three formulas we obtain: 



1> (a, b; e^z) = ^ (a, b; z) e~ 2mb + cf> (a, b; z) - - - _ - e-" b . (52) 
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Set 



rft+t) ^(A + i.7 + i;C)^ 



-C/2 e i7r 7 /4 



12 



-C /2 ^ (a + |, 7 + 1;C) e-^V^/ 4 . 



Then from (51) and (52) if follows that for £ G Ts, 



11 



(0 = (e 2m cf 2 4> (A + 1, 7 + 1; e 2 -c) e^V^^gfl 
= e^(W) (C), 



and 



12 



(0 



e -C/2 e -i7r 7 /4 



27r^e-"e-' ,r7 / 4 r(7+l) 

r(A-^)r( 7 +i)r(i-A+^)e i ' r T/ 4 



_ e -i7rA _|_ e «7rA e -i7r7 



*i0 (C) + e" (7 ~ 27 ~ 2) (*l2) (C) 



(0 



*n (0 + e- W7 (*i2 



in accordance with (49). Analogously, we can satisfy the second row of (49) if we take 



21 



22 



r(i + A + 

r( 7 + i; 

r(i + A + j) 
r(|-A) 



^r /2 0(l + A+^ 7 + l;C) 



C 7/ V(l + A + |,7 + l;C) e-^e-~i. 



By construction, \I/ satisfies the jumps relations in (^2). Using formulas (7.26), (7.27) and 
(7.29) from |13[ appendix], we can write explicitly the matrix \I/ in all regions. Since the local 
behavior of ip(a, b; z) depends only on the value of the parameter b, by construction, all rows 
of \I/ have the same asymptotics as £ — > 0. Hence, it is sufficient to analyze the first row. 

From formulas (13.5.5) and (13.5.12) from jT] if follows that for £ G ®, \I/ has the behavior 
described in (^3), as £ — > 0. Indeed, for 7 > 0, 



11 



for 7 = 0, 

and for — 1 < 7 < 0, 



*n = 0(l), *i2 = 0(lnC); 
*n = 0(c /2 ), *i2 = 0fC 7/2 V 
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Analogously we can check that \l/ satisfies (^3) in all regions of the plane. 
Finally, using formula (13.1.22) from [T], 



<f>(a,b;C) tp(a,b;C) 
0'(a,b;() i>'{a,b;0 



r (b) e< 



C b r(a)' 

as well as the differential relations (13.4.23) and (13.4.10) from [I], we easily get that 



</ 2 e^ -C^V(a,6;C)e-^ 2 e-T 
'WC^Ha + 1, b; C) e-^eT r /ffi?> C*^ (« + 1, * e^/ 2 



^C^0(a,6;C)e" 



e 4 



1. 



r(6) s vvtj.,!/^^ r(-(i+a-b)) 1 
This implies tha det \& = 1. and, by construction, det = 1, which concludes the proof. 



In order to construct the analytic function E n in (29) we need to study also the asymptotic 
behavior of \I/ at infinity. Let us introduce the notation 

(*+i)„ (*-§)„ 



v n = v n (A) 



ni 



-A+l) 



r(|-A) 



H-A)r(i + A)" r(i + A + i)' 



observe that 



t-a = Tx> (-A) = d„ (A), and v\ = A 



(53) 
(54) 

(55) 



Lemma 4. As ( -> oo ; ( G C\r, 



1 



(-l) n w n nr A fr. 



+ o ICI 



C-^e-^M" 1 ^) (56) 



with v n defined by (53), t\ defined by (54), A = zlog(c)/7r, and 



M(C) 

where we use the main branch of 



e 4 

P.4 



?7ri<T3 A7ri<J3 



-1 

o i N 
-1 o. 



I < argC < 7r, 
7T < argC < | [ , 

-f < argC < 0, 



< argC < | , 
e -AlogC y^/j i/jg Clt £ along il 
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Proof. We use the classical formulas (13.5.1) and (13.5.2) from [T] for the confluent hyper- 
geometric functions. If we take 6 = 7 + 1, and multiply <f> and if) by (z 7 / 2 e _2//2 ), using (38), 
we have that, as \z\ — > 00, 



G(a,r,z) 



r(7+i) ~a- 7 /2 

r(a) z 

r (7+l) ^aTTt r7 /2-a 

r( 7 +i-a) e Z 



JJ-1 ( 7 +l-a) n (l-o) r , 



I _u V^-R-l («)n(a-7)» , a) / I . 
1 -T Z^n=l (-l)"n!«" +U \ \ Z 



-R 



-R 



H(a,r,z) = z 



7/2-0 



e^ 2 , Rez>0, 

e~ z / 2 , Rez<0, 
(57) 



1 , V" M"^-^ I g f U| 



n=l 



-z/2 



(58) 



Replacing these expansions in the expression for \& for £ 6 ®, | < arg C < x an< ^ < - 



arg 



(e-^C) < we get for |C| - 00, 



*(C) 



C" 



1 + E*=i 



1 (*+2)„(*-2)„ 



(-l) n n!C 



+0 ici 



-r(i+A+g) 1-A 
V r (5-^) ^ 



e -C/2 e A7ri 



r(i-A+i; 

r(i+A) 



l 



-l+A 



1 , ( 1 +A+I), 1 (1+A-|) 
1 Z^n=l (-l) n n!C n 

-1 (l-A+i) n (l-A-|) T 



+0 ici 



1 + E^=i 



1 + E*=i 



(-l) n (-l)"C™n! 



+ Id 



1--R 



(-l) n (-l) n C n n! 



" +0 id 



=C/2 



e -C/2gA7ri 



=C/2 



e 4 



0-3 



which can be rewritten using notation (53)-(54) as 



1 + 



En=l ("I)" P 
sr^R— 1 / -i \n m 



R—l v n 



n=l 



+ id 



This yields (56) for 7r/2 < £ < 3tt/4; this expansion is also valid for £ € ©. A comparison 
of (42) with (43) shows that the behavior for ( G (D, 7r < arg£ < can be obtained from 
the expansion in © by multiplying by e*2 7r(73 , which again yields (56) for 7r < ( < 5ir/4. It is 
easy to see that asymptotics in © is also valid in ®. 

Using (45), (57), (58) and comparing the expression for \I/ in ® and ©, we conclude that 
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for Ce ®, -f < argC < -f (f < arg (C) e ni < & and ReC > 0), as |C| oo, 



*(C) 



/ r(i-A+j) 
r(i+A) 



-l+A 



1 + V^ 1 ( 1 ~ A+ f)n( 1 ~ A ~i)n i fo/V|-R 

1 ~r /Lm=l (-l)"(-l)"C n n! " l " Ly U t 'l 



1 + £„=l 



1 (I-A) n (-|-A) n 



(-l)"(-l) n C™n! 



+o ici 



g-ATTigC/2 



-r(l+AH 



1 I V^ 1 ( A+ i)n( A ~j) n , ^ / 1 ^-R 

2—jti = \ (-l) n n!C n "r^ll^l 



-C/2 



r (i" A ) 



-l-A 



1 T Z^n=l (-l) n n!C n "T^U^I 



-1 

1 

-1 

1 

R-l 



R-l 



n=l s 



m\v n 



-{-l) n v n 
{-l) n nT\v r . 



"C/2 



e * 



^3 



+ c Id 



-1 

1 



1 



(-l) n v n nr x vr, 



+ Id 



-ii 



\ 1 ! C 

£— A(J3 e ^7ri<T3 e — J<X3 



-1 

1 



This expression is valid in © as well. Finally, comparing (46) with (47) we see that the 
behavior for £ G ©, < arg£ < ^, corresponds to that in © times the constant factor e*2 7rfT3 , 
which yields (56). Since the asymptotics for £ G ® is the same than in ®, this concludes the 
proof of Lemma. 

Now we are ready to build in (29). Using the properties of ip we define an analytic 
function / in a neighborhood of Xq, 

de f I 2i arccos xq — 2 log <p (z) , for Im z > 0, 



/(*) 



2i arccos xq + 2 log tp (z) , for Im z < 0, 



(59) 



where we take the main branch of the logarithm. Using that ip + (x) <p- (x) = 1 on (—1, 1) 
we conclude that / can be extended to a holomorphic function in C \ ((— oo, —1] U [1, +oo)). 
For \z\ < 1 we have 



/(*) 



2i 



(Z - X ) + O ((z - Xq) 2 ) , as Z 



Xq. 



Hence, for 5 > sufficiently small, / is a conformal mapping of U XQ . Moreover, since 
then 



p + (x) = x + iVl-x 2 = e iaT ' 
f(x) = 2i (arccos xq — arccos x) , x G (—1, 1), 



(60) 

(61) 
(62) 
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so that / maps the real interval (— 1,Xq) one-to-one onto the purely imaginary interval 
(2i(arccos xq — ir), 0), as well as (xq, 1) one-to-one onto the purely imaginary interval (0, 2i arccos xq). 
We can always deform our contours Y,^ close to z = xq in such a way that 

f(E x nu X0 ) cr 4 , f(z 2 nu xo ) cr 6 , /(E 3 nC4 ) c r 2 , /(s 4 nf/ I0 )cr 8 
f(E 5 nu X0 ) cr 3 , /(E 6 n£/ ao ) cry 

With this convention, set 



and, we define 



( = nf(z), zeU X0 , 



(63) 



(64) 

By (\]/l)-(\]/3) and (60), this matrix-valued function has the jumps and the local behavior at 
z = xq specified in (32)-(35). Taking into account the definition (59) we get that 

e nf(z) = ^( Xo ) ^ 2n ( z ), for ilmoO, 



and for [nf (z)] we take the cut along (— oo,xq]. Since 



[/(z)] A = |/(z)| A e X p 



log c 



7T 



arg(/(^)) 



straightforward computations show that 

'l/OzOlV 1 / 2 , forx <x<l, 
\f{x)\ X c- 1 ' 2 ^ 1 , for - 1< x < x Q , 

where we assume the natural orientation of the interval. 
In order to satisfy (Po3) above, we define 



(65) 



B n (z) = N(z)W(z) 



o r 
-i o. 



(nf(z))^^(x )e 



(nf(z))^^(x ) 



o r 
-i o. 



(nf(z)) x ^^(x )e- 



if z € Qf, 
if z G Q R , 

if z e Q^, 

if ze o 1 :. 



(66) 



By construction, E n is analytic in C4,. \ (K U S5 U Eg). Furthermore, by (N2) and (28), for 

X £ (Xq- 5, Xq) U (x , Xq + 5), 

W- { x y a[i n: 1 (x)n + (x) w + (x) CT3 







W-(x)W+(x) 
W-(x)W+(x) Q 
Wc, 7 (z) 

for =b X > Xq ; 



C 

-c^ 1 



20 



and, by (26 ), for z G Eg H {7 X0 (oriented from above to bellow) and for z € E5 n (oriented 
from bellow to above) we have, 



(z)-* 3 nz 1 (z)n+(z)w + (z) c 



W + {z)/W„{z) 









W-(z)/W+(z] 



e 2 



From ( 65 ) and ( 66 ) it follows that 

E~i (z) E„ + (z) = I, for z G ?7 X0 \{x } . 

In this form, xo is the only possible isolated singularity of E n in U xo . The following proposition 
shows that this is in fact a removable singularity of E n : 



Proposition 5. 



lim E n (z) 

z—>x 



V2 



3 — i arcsin(iEo)/2 i arcsin(:ro)/2 



2 00 I gi arcsin(xo)/2 g— j arcsin(xo)/2 



wi/i 7? n defined by 



def 



log C 



7T 



log y^nyj 1 — Xq^ +narccos(xo) 



77T 



$(x ) 



(67) 



and $ given by (22). /n particular, E n is analytic in U XQ . 



Proof. Since E ra is analytic in a neighborhood of xo, it is sufficient to analyze its limit as 
z — > xq from the first quarter of the plane, z £ Q5- By (24) and (60), 



lim D(z,H c )/(zr A = lim c 1+ ^ lo s( z / 2 )-^ M/tt) = c 3/2 ( 4 /1 _ x 2 



On the other hand, by (21 ) and d27j) (notice that u>i i7 defined in ( 12 ) coincides with w defined 



in (11)) 



lim D(z,w) W(z)- 1 = c -l/2 e i*(^) e i7^/2_ 

z — x 



Summarizing, 



lim D(z,w c ~,)~ x W (z) f (z) X 



4^1 



-i$(xo)— iyir/2 
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By (14b and (|66j), if z G (Imz > 0) 



with 



E n (z) = L»SA(z)m n (z) CT3 , 



D(z,w cn ) 



(68) 
(69) 



where ?? n is defined in (67). 

Gathering the limits computed above and using that 



lim A n (z) = e - iarcsin (^)/2 = lim A , z } 



if? 



and by definition of rj n , the statement follows. 



Therefore, by construction the matrix- valued function given by (29) satisfies condi- 
tions (Pol)-(Po4). Moreover, it is easy to check that 

det P XQ (z) = 1 for every z G U XQ \E. 

At this point all the ingredients are ready to define the final transformation. We take 



R(z) 



ef f S (z) N" 1 (z) , z G C\{S U [/_! U ^ U *7i} ; 



[S^P-^z), z eUj\ E, j G {-l,x ,l}. 
R is analytic in the complement to the contours Er depicted in Fig. [4] 



(70) 




Figure 4: Contours "Sr. 



Let 



d = f E R \ (9t/_! U dU X0 U at7i) 
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R satisfies the jump relation 



R+(z) = R_(z)(I + A(z)), z £ E fl , 



with 



N( S ) 



u> Cj7 (s) 1 <p(s)~ 



2» 



N- 1 (s) - 1, for s G S 



out. 
R ' 



l P C ( s ) ( s ) _I > for S GaC/ c , ;€{-l,x ,l}. 

Standard arguments show that A has an asymptotic expansion in powers of 1 jn of the form 



A (s) ~ — j as n — > oo, uniformly for s £ Hr, 



k=i 



(71) 



and, for k G N, 



A fe (s) = 0, forsGS^*. (72) 
So, it remains to determine A/% on dU Xo - Here we find explicitly only the first term, Ai. 



Using (|14|, fl25J, p9\, (|59|), (|56|), (|63|, (|66| and (|55j), we obtain 

1 



A(s) = E n (s) 
Let us define 



(A 2 "7 2 /4)/-l r A 



n/(s) 



-n i 



+ 



E n (s), s £ dU Xo , n — > oo. 



A , , def (A 2 - 7 2 /4) w , > 



Using that by ( 66 ) , 



E n (a) = P(a) l^ + (x Q ) n n x ) " = F(a) (e 



3 inarccos(:ro)g^: logn 



(73) 



with 



F/ s. del 



'N (a) W (s) CT3 c ff3 e ±2 ? ff3 / (s) A,T3 
N(s)T^(s 



CT3 ( °i o) e±2?,73 ^( s ) A<T3 ' if Ims < 



if Im s > xo , 



where we take db for d=Res > xq, we conclude that, for s £ dU X0 , Ai(z, n) is uniformly 
bounded in n; indeed, F does not depend on n and 



% log n 

in arccos xn „^r l — 
e C 27T 



1, Vn G N. 
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So Ai in (73) is genuinely the first coefficient in the expansion (71). 

Similar analysis can be performed for A^ (-,n), k > 2, taking higher order terms in the 
expansion of ^ in ( 56 ) . 



The explicit expression (56 ) and the local behavior of / show that Ai (s, n) has an analytic 
continuation to U xo except for xq, where it has a simple pole. Again, similar conclusion is 
valid for other A^ (s,n), except that now the pole has order k. 



Like in p2 Theorem 7.10] we obtain from (71) that 

R (z) ~ I + ~. ~ — , as n — :\ . 



(74) 



uniformly for z £ C\ {dU-i U dU XQ U dUi}, where each Rw) (z) is analytic, uniformly bounded 
in n, and 

R (3) (z,n) = oQ as ^ oo. 

Since RW is analytic in the complement of dU-i U dU Xo U dU\ (see (72)) and vanishes at 
infinity, by Sokhotskii-Plemelj formulas, 



R (1) (z,n) = -^ / 

2m JdU-T. 



Ai(s,n) 



ds. 



udU XQ \JdU! s ~ z 



Recall that Ai can be extended analytically inside t/j's with simple poles at ±1 and xoi let us 
denote by (n), B^ (n) and (n) the residues of Ai(-, n) at 1, —1 and xo, respectively. 
Then the residue calculus gives 



R (1) (z,n) = < 



( (n) BW (n) C^(n) 

Z — 1 Z + 1 Z — Xq 

AW ( n ) BW (n) C« (n) 



z- 1 



+ 



Z+l 



+ 



for z EC\ {U-i U U XQ U Z7i} ; 
Ai (z, n) , for z G U C4 U Ui. 

(75) 

Residues A^ 1 ) (n) and B^ (n) are in fact independent of n; they have been determined in 
[13 Section 8]: 



z- XQ 



(n) = A& 



4a 2 - 1 



D a 3 

16 00 



-1 t 
i 1 



^OO 1 



B a. (n ) = B a» = ^!_i DS ^ 



(76) 



0":i 



(notice however that the constant is different with respect to [13] )• The value of the 
remaining residue (n) is given in the following 
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Proposition 6. If we denote 

C« (r 

then the entries are given explicitly by: 



n 



C{\\n) C&\n) 



^(1)/ \ /log 2 c 7 2 \ /log 2 c 7 2 . 



0$ (n) = il& [ ^ + £ - J ^ + £ cos (arcsin (x ) - n ) j ( 78) 



<4? (") = 7^ ( ^ + t + \/^ + t «» (—in (xo) + *„) ) (79) 



1)2, \ 2vr 2 8 V 4vr 2 16 



where 



9 n = 2% + (81) 



mi/i ?y n defined by (67) and ? = — 2 argT + A) — arg {% + A) . 



Proof. Taking into account (60) and ( |73[ ) we conclude that the residue (n) of Ai (z,n) 
at z = xo is given by 

cm („, = (^-7V4)v^ En (n) ( -i_ ? ) E - t (lo) . (82) 

Since E n is analytic in a neighborhood of x$ (see Proposition]^]), 

1 / — i arcsin(xo)/2 i arcsin(a;o)/2 

so that 

^ / g— j arcsin(a;o)/2 gi arcsin(a;o)/2 \ 

E n (xq) = — ^=— -j===e %T]n(JZ ^ e iarcsin(a;o)/2 g -i arcsin(x )/2 J Axf 3 - (84) 



E n (x ) = hm E n (z) = -—==^D£ / _ ejarcsin(xQ)/2 e _ iarcsin(3;o)/2 ) e™, (83) 
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From (83) we obtain 



n 



(A 2 - 7 2 /4) 



D a ' A 



-i arcsin(a;o)/2 i arcsin(a;o)/2 



i arcsin(a'o)/2 



— e 

"I T A 

-rx 1 



e 

— irinCTi 



i arcsin(a:o)/2 

g— j arcsin(xo)/2 
gi arcsin(a;o)/2 



i arcsin(xo)/2 

1 D~ a3 

— tarcsin(x'o)/2 ; oo 



Using formulas (6.1.28), (6.1.23) and (4.3.2) from [T] and (54) we can rewrite 



rQ + A) 



;i + A)r(2 + A ) 



r(i + A) 



+ A)r(2 + A ) 



^ T V4 + |A| ; 



where 



? = arg 



r(i + A) 



(i + A)r(i + A) r 



Then, 



C (D ( 



(A 2 -7 2 /4) 



x - 
-i — i 



Vt 2 /4+|A| 2 
cos(arcsin(a;o)+Sn) 

V7 2 /4+|A| 2 



cos(arcsin(a;o)— 6 n ) ' 
V7 2 /4+|A| 2 

Sin(-fln) 



D' a:i 



V7 2 /4+|A| 2 



v2 _ 



We can simplify this expression using that ^A' j 
and this settles the proof. 



log 2 c _i_ t 2 

^■2 "T A 



3. Proof of Theorem 1 

Unraveling the transformations Y->T-»S^Rwe can obtain an expression for 
Y. Repeating the arguments in [8] (see also [H], [HJ Section 3] and [T7]) we see that the 
recurrence coefficients ^ are given by 



lim 

z— >oo 



£> 2 A / l 

" " + zR 12 (z, n) zR 2 i (z, n) + 



2% 



lim z (1 - Rn (z, n + 1) R 2 2 (-2, ra)) 



(85) 
(86) 



Taking into account the expression for R^ in (75), as well as (76) and Proposition [HJ we 
obtain for a n : 



1 1 / 7 2 log 2 c 
4 ~ nV 16 + 4vr 2 



cos (arcsinxo) cos (0 n ) + O 



n — * oo, 
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where 9 n is given by (81 ). It also can be rewritten in the form 

2 log c 



7T 



log ( An\/1 — Xq ) + 2n arccos xq — 6, 



with defined in Q. This proves Q. 
Analogously, 



2 2 

+ 16 ( sin ~ sin n ) / i 



By (81), 



b n = ILj : L + o(JL.) n _> oo. 

n \n z 



log c 

Wi - = 2 arccos(xo) + 2 log ( 1 + - ) , 

7T V Tl 



and 

1 



sin ^n+i — sin 9 n = sin (Q n + 2 arccos xq) — sin 6 n + O 



n 



2cos( 9 „ + 8I ccosx )sin(arcco S , 0)+O |I], ,w =c. 



Thus we obtain 



1 /log 2 c 7 2 ( \ 

—\\ — 5 1 — — sin (arccos xq) cos (6 n + arccos xq) + O \ —r 

n V tt z 4 \ti^ 



which proves ([5]). 
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